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Abstract 

An acoustic prediction capability for supersonic 
axisymmetric jets was developed on the basis of 
OVERFLOW Navier-Stokes CFD (Computational Fluid 
Dynamics) code of NASA Langley Research Center. 
Reynolds-averaged turbulent stresses in the flow field are 
modeled with the aid of Spalart-Allmaras one-equation 
turbulence model. Appropriate acoustic and outflow 
boundary conditions were implemented to compute time- 
dependent acoustic pressure in the nonlinear source-field. 
Based on the specification of acoustic pressure, its 
temporal and normal derivatives on the Kirchhoff surface, 
the near-field and the far-field sound pressure levels are 
computed via Kirchhoff surface integral, with the 
Kirchhoff surface chosen to enclose the nonlinear sound 
source region described by the CFD code. The methods 
are validated by a comparison of the predictions of sound 
pressure levels with the available data (Seiner et al., 1992) 
for an axisymmetric turbulent supersonic (Mach 2) 
perfectly expanded jet. 

Nomenclature 

c =sound speed, ^yg c R u T IW , ft/s 
dj = nozzle exit diameter, ft 

F = total thrust, lbf 
/ = frequency, Hz 

g c = gravitational constant, (lbm-ft/lbf.s 2 ) 

/ = specific impulse (sec) 

M = Mach number, Vic 
m = instability mode 
p = static pressure, lbf/ft 2 

p = acoustic pressure disturbance (instantaneous 
pressure), p - p m , lbf/ft 2 
Q = vector of disturbance variables 

Re = Reynolds number, pujdj / p 

R u = universal gas constant, lbf.ft/(lbm-mole R) 

r = radial distance 

rj = nozzle exit radius, ft 

St = Strouhal number, fdj luj 
T = temperature, R 


/y a 

V = velocity, m/s 

u = axial velocity, ft/s 

u c = jet center-line velocity at exit, ft/s 

w= radial velocity, ft/s 

W = molecular weight, lbm/lbm-mole 

x = axial distance from the nozzle exit, ft 

z = radial distance from jet axis, ft 

Greek Symbols 

fi = dynamic viscosity, lbm.ft/s 
p = density, lbm/ft 3 

Y = isentropic exponent 

co = circular frequency, radians/s 
y/ = angle from jet axis, deg 

Subscripts 

c = centerline 
e = nozzle exit 
7= jet 

t = stagnation or total, or nozzle throat 
oo = ambient 

Superscripts 

prime : disturbance from the mean value 
1. Introduction 


Acoustic loads in a launch vehicle environment 
represent a principal source for inducing structural 
vibration and may be critical to the proper functioning of 
vehicle components and ground support structures and 
equipment in the immediate vicinity of the launch pad. A 
knowledge of acoustic loads, including the overall sound 
pressure level (OASPL), sound pressure level (SPL) 
spectrum, and the distribution (or correlation) of surface 
acoustic loads, is necessary to provide the input for 
vibroacoustic analysis and evaluation of structural 
integrity. In the design of launch vehicles, it is highly 
desirable that data on acoustic loads (near-field and far- 
field noise levels) be generated both analytically and from 
testing of small-scale and full-scale models. Since full- 
scale acoustic and vibration testing is often cost 
prohibitive, the option of small scale testing combined 
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with analysis methods remains as a practical alternative. 
Accurate characterization of acoustic loads on launch pad 
structures thus proves to be a formidable challenge 

Noise from subsonic jets is mainly due to turbulent 
mixing, comprising the contributions of large-scale and 
fine-scale structures (Lighthill 1,2 ). The turbulent mixing 
noise is mainly broadband. In perfectly expanded 
supersonic jets (nozzle exit plane pressure equals the 
ambient pressure), the large-scale mixing noise manifests 
itself primarily as Mach wave radiation, caused by the 
supersonic convection of turbulent eddies with respect to 
the ambient fluid. In imperfectly expanded supersonic 
jets, additional noise is generated on account of 
broadband shock noise and screech tones. 

The prediction of noise generation and propagation 
from exhaust jets by computational methods entails 
certain numerical requirements (Tam 3 ). Although the 
existing CFD codes can predict the noise generation in the 
nonlinear source-field, their application to acoustic 
predictions for obtaining near-field and far-field sound 
levels results in inaccuracies due to numerical dispersion 
and dissipation over long propagation distances. 
Specialized Computational Aeroacoustics (CAA) methods 
were developed to resolve these numerical issues 
satisfactorily (Tam and Webb 4 , Lele 5 ). These 
computational methods for direct calculation of noise in 
the near- and the far-field, however, require prohibitively 
large computational resources. 

This situation motivated in recent years the 
application of CFD codes to model the source region 
(containing the sound sources), with the near- and far- 
field sound predicted by either Lighthill’s acoustic 
analogy (Lighthill 1,2 ) or the Kirchhoff surface integral 
(Kirchhoff6). In Lighthill’s acoustic analogy the Navier- 
Stokes equations are rearranged in wave-like form with 
Lighthill's stress tensor acting as the source term, and the 
near-field and far-field sound is obtained via a volume 
integral. This volume integral again requires large 
computational resources. On the other hand, the Kirchhoff 
method requires only a surface integral, with the 
Kirchhoff surface enclosing the nonlinear source region, 
and is chosen in a region where the linear wave equation 
is valid. Additionally, the Kirchhoff method does not 
suffer dissipation and dispersion errors when the near- 
field and far-field sounds are directly calculated with the 
CFD codes. By restricting the use of CFD methods to the 
nonlinear near-field region for source identification and 
employing Kirchhoff integration for the linear region, 
computational requirements are greatly reduced. 

The ability to accurately predict the turbulent mixing 
noise sources remains as a key requirement in the overall 
acoustic analysis. Both the large-scale turbulence 
structures (instability waves) and the fine-scale turbulence 
represent the mixing noise sources. The large turbulence 
structures are responsible for Mach wave emission that is 


directional, and the fine-scale turbulence mixing noise is 
fairly uniform independent of direction. In the past the 
two-equation k — £ turbulence model is generally used 
to model the sound sources. However, this model entails 
large computational requirements (two additional 
equations) and is often plagued with convergence 
problems in flows with shock waves and strong gradients. 
Furthermore, the universal constants in the turbulence 
model are often altered to fit the theory with the acoustic 
measurements. On the other hand, the simple algebraic 
mixing length type turbulence model cannot accurately 
describe the growth of the mixing layers. Thus a one- 
equation turbulence model seems to be an attractive 
choice for engineering applications. This model generally 
results in faster convergence and requires less 
computational resources. 

The objective of this investigation is to predict 
supersonic jet acoustics on the basis of CFD/ Kirchhoff 
formulation in conjunction with a one-equation turbulence 
model. This report summarizes the development of a 
computational procedure in which the OVERFLOW CFD 
Navier-Stokes code (Buning et al. 7 ) predicts the nonlinear 
source-field within the jet and application of Kirchhoff 
surface integral for acoustic radiation. 

2. Analysis 

This method of analysis is based on the application of 
a CFD code for identifying the noise sources in the 
source-field and Kirchhoff surface integral for the 
propagation of sound radiation to the near-field and the 
far-field. 

2.1 CFD Analysis 

2.1.1 Numerical Algorithm 

For the CFD analysis, the OVERFLOW code 
(Buning et al., 1998) of NASA Langley is used for 
computing the instantaneous flowfield. This code solves 
the three-dimensional compressible turbulent flow 
Navier-Stokes equations in generalized coordinates. The 
one-equation turbulence model of Spalart-Allmaras 8 is 
considered. 

2.1.2 Boundary Conditions 

To simulate the exterior (infinite boundary) acoustic 
problem, we need to prescribe outflow and radiation 
boundary conditions at the boundaries of a finite 
computational domain. These boundary conditions must 
be sufficiently transparent to the outgoing disturbances so 
that they exit the computational domain without 
significant (nonphysical) reflections. It is assumed that 
the boundaries are sufficiently far from the sources. 

The OVERFLOW code, in its present form, is 
designed primarily for the prediction of steady or 
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unsteady turbulent flowfields and used widely to model 
aerodynamic flows. It does not however provide the 
necessary time-dependent boundary conditions for 
handling the reflection-free acoustic propagation. 
Therefore appropriate modifications to the code were 
made to provide a set of time-dependent reflection-free 
boundary conditions (includes periodic inflow, outflow, 
and radiation). 

a. Inflow Boundary Conditions . 


2.2 Kirchhoff Code 

For the acoustic radiation, the Kirchhoff code 
YORICK developed by Pilon and Lyrintzis 11 of Purdue 
University was considered. The Kirchhoff surface, 
enclosing the nonlinear source region, is chosen in a 
region where the linear wave equation is valid. Outside 
the Kirchhoff surface, the flow is governed by the linear 
equation 


At the supersonic inflow, all data are specified, such 
that all time variations of the characteristic variables are 
set to zero because all waves are incoming waves. 




v V--Thr + “~— p'f 0 


dt 


dx. 


(2a) 


b. Outflow Boundary Conditions . 

In general at the outflow boundary, the outgoing 
disturbances consist of acoustic, entropy, and vorticity 
waves. Here we follow Thompson’s approach 
(Thompson 9 ) based on one-dimensional characteristics 
analysis (Riemann variables). The amplitude of the 
outward propagating waves at the boundary are specified 
entirely from the solution within the computational 
domain, while the amplitude of the inward propagating 
waves is prescribed at the boundary. For details on the 
method, the original reference should be consulted. 

c. Radiation Boundary Condition . 

At boundaries where there exist only outgoing 
acoustic waves, Tam and Webb 4 developed a set of 
radiation boundary conditions based on an asymptotic 
analysis (valid for large distances). These boundary 
conditions, originally formulated in a spherical coordinate 
system, when applied to a cylindrical coordinate system, 
are given by (Hixon et al. 10 ) 


where u ^ is the freestream (ambient) velocity. In the 
special case of stationary ambient ( = 0 ), it reduces to 
the simple wave equation. A solution to the pressure field 
can be expressed by the surface integrals as (Morino and 
Tseng 12 , Lyrintzis & Mankbadi 13 ) 

p'( X ,y, Z ,t) = --L 

4 it 

" P_ d'o ! 1 dp' | 
r 0 2 dn 0 r 0 dn 0 

1 dp' f dr 0 M fop ^ 
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In the above equations, Q' is the vector of disturbance 
variables, x the axial distance from the origin, and r the 
radial distance from the jet axis. The subscripts on 
Q' indicate partial derivatives. 


In the above equations, is the freestream Mach 

number, c w the speed of sound in the freestream, the 
prime a point on the Kirchhoff surface, T the retarded 
time T = t -t\ n 0 the outward normal to the Kirchhoff 
surface S 0 . Thus the pressure at any instant in the region 
outside the Kirchhoff surface can be expressed in terms of 
the information prescribed on the Kirchhoff surface. The 
required data on pressure, and its normal and temporal 
derivatives on the Kirchhoff surface are taken from the 
CFD solution. 

3. Validation Test Case 


The CFD/Kirchhoff methodology formulated above 
is tested by comparing the predictions with the well- 
known acoustic test data of Seiner et al. 14 for Mach 2 
perfectly expanded supersonic jets of air. The nozzle exit 
radius is 0.04572 m. Although the data cover a wide 
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range of jet total temperature (using propane burners), 
detailed validation is performed here for a jet total 
temperature of 755 kelvin (K). Validation at other jet 
total temperatures will be reported in a later report. Table 
1 shows the data at the nozzle exit. The jet exit velocity is 
822 meters per second (m/s), and exhausts into quiescent 
surroundings at a temperature of 280 K. Figure 1 shows 
the location of the far-field measurements of the OASPL 
with which comparisons with the CFD results are 
reported. 

4. Numerical solution 

4.1 CFD Solution 

4.1.1 Axisvmmetric Grid 

For the CFD computations, an axisymmetric grid of 
size 200 x 100 x 3 (200 grid points in the axial direction, 
100 grid points in the radial direction, and 3 in the 
azimuthal plane) is considered (figure 2). A grid length 
of 60 jet radii and a grid radius of 10 radii are considered. 
The grid is clustered in the radial direction near the nozzle 
wall to resolve the shear layer. Because of the perfectly 
expanded condition, uniform grid is considered adequate 
in the axial direction. 


about 600 iterations for convergence) has served as the 
starting solution for the time-accurate run. The solution is 
obtained on an IRIX workstation (SGI Indigo machine). 
About 15 hours of computing time are required for the 
solution to achieve a periodic state. 

4.2 Acoustic Solution 

After a periodic state has been established, the 
appropriate data from CFD, corresponding to a cycle, is 
communicated to the Kirchhoff code. The Kirchhoff 
surface is a cylindrical surface coinciding with a grid line. 
The radius of the Kirchhoff surface is taken as about 6 
radii from the jet axis. Only the lateral surface of the 
cylinder is taken into account, and the cylinder ends are 
neglected due to the effects of nonlinearities. The data to 
be specified on the Kirchhoff surface include the 
instantaneous pressure, the pressure-time derivative, and 
the pressure gradient normal to the surface. The 
Kirchhoff code then computes the OASPL and the 
pressure-time signals in the near- and the far-field. From 
the pressure-time signals, it is possible to compute the 
SPL spectrum at any location with the aid of the Fast 
Fourier Transform (FFT). 

5. Results and Discussion 


4.1.2 Flowfield Solution 


5.1 Jet Mean Flowfield 


At the inlet, the flow variables are specified based on 
the nozzle exit conditions. In addition, a periodic 
disturbance (axisymmetric mode) of a single frequency is 
also prescribed as follows: 


u'(z)= Asin(cot)cxp 


-ln(2 


z-h 


( 3 ) 


where CQ is the circular frequency of the disturbance and 
A the amplitude of disturbance. Based on the work of 
Mankbadi et al. 15 , values of 6 = 0.1 and h = 0.78361 are 
considered. The amplitude A is taken to be of the order 
of 0.01. A value of Sr = 0.2 (based on typical peak SPL 
frequency for supersonic jets) is considered here for 
obtaining the disturbance frequency CD . 

Appropriate time-dependent boundary conditions 
were applied to ensure reflection-free boundaries. An 
outflow boundary condition of Thompson type 
(Thompson 1990) is applied at the outflow boundary, 
which maintains the mean static pressure at the ambient 
value. An acoustic radiation boundary condition of Tam 
and Webb 4 is applied at the lateral boundary. 

The CFD solution converges after about 8,000 time- 
step iterations (the code is run in a time-accurate manner) 
before a periodic state is established. An initial steady- 
state solution (based on local time-stepping, requiring 


Figure 3 shows the computed mean velocity contours 
in the jet. The potential core and the mixing of the 
turbulent shear layer are clearly displayed. A comparison 
of the predicted centerline mean axial velocity with the 
data is presented in figure 4. Both the data and the CFD 
solution suggest that the jet potential core extends to 
about 15 jet radii from the nozzle exit plane. The jet 
potential core may be defined as the axial distance where 
the jet centerline velocity remains at about 99% of the jet 
exhaust velocity. It is known that the major noise 
production region occurs near the end of the potential 
core, where the most highly amplified instability wave 
reaches its maximum growth (Seiner et al. 14 ) 

In the velocity decay region past the core, the CFD 
solution predicts a slower growth relative to the 
measurements. This discrepancy is attributable to the 
accuracy of the Spalart-Allmaras one-equation turbulence 
model. Also, the effects of compressibility are not 
appropriately taken into account in this turbulence model, 
as it was originally developed for incompressible flows. 
Data on compressible shear layers indicate that 
compressibility tends to reduce the growth rate of the 
shear layers (Goebel and Dutton 16 ) relative to the 
incompressible case, with an increase in the jet convective 
Mach number. For example, before the potential core 
closes, a jet convective Mach number M c , which may be 
defined by (Freund et al. 17 ) 
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The reduced transverse turbulence intensity is believed to 
be responsible for the observed reduction in the spreading 
rate for compressible mixing layer. By taking into 
account the effects of compressibility on the turbulent 
kinetic energy and its specific dissipation rate into a two- 
equation turbulence model, the reduction in the mixing 
length growth rate is found to agree with the data 
(Kandula and Wilcox 18 ). 

Figures 5 and 6 respectively display the mean static 
pressure and mean Mach number along the centerline of 
the jet. The pressure oscillation along the centerline is 
relatively small and shows the accuracy of the boundary 
conditions. The centerline Mach number variation 
follows the trend of the centerline velocity, as is to be 
expected. 

Calculations show that with an increase in jet 
temperature, the potential core length decreases. These 
results are in qualitative agreement with the data of Seiner 
et al. 14 . As the jet temperature increases, the convective 
Mach number increases, and the Mach wave radiation 
pattern is modified, with the OASPL values generally 
increasing with an increase in jet temperature. 

5.2 Jet Instantaneous Pressure and OASPL Contours 


Figure 7 displays the instantaneous pressure contours 
in the source-field, showing the directivity of Mach wave 
emission. The lack of significant reflection near the 
boundaries indicates that the outflow and the radiation 
boundary conditions implemented in the CFD code are 
satisfactory. A display of the OASPL within the jet, 
characterizing the nonlinear source-field, as directly 
computed by the CFD code, is presented in figure 8. 
Inside the jet shear layer, the pressure levels are high (as 
much as 180 dB), with the highest OASPL occurring 
about 7 jet radii downstream of the exit. 

5.3 Near- and-Farfield Acoustic Solutions 


Figure 9 presents the OASPL contours in the near- 
and the far-field as computed by the CFD/Kirchhoff 
formulation. It is seen that the OASPL ranges from 80 dB 
to 130 dB in this region. The directivity of the Mach 
wave radiation (emission) is evident. It is generally 
accepted that the mechanism for Mach wave radiation is 
strongly connected with the amplification of instability 
modes, and this amplification scales with the jet 
convective Mach number. 


measurements of Seiner et al. 14. The predicted peak angle 
of emission of about 128 degrees is in good agreement 
with the data, which indicates a peak at 127 degrees. 
Also, the peak value of the OASPL is about 7 dB less 
than the measured value. There is an underprediction of 
the OASPL at the off-peak locations both upstream and 
downstream of the peak. This underprediction is 
primarily due to the axisymmetric simulations 
(axisymmetric disturbance mode m = 0 instability) of the 
present investigation. It is impossible to develop a three- 
dimensional helical mode instability with axisymmetric 
simulations. Measurements show that near the nozzle 
exit, axisymmetric modes seem to be dominant, and 
farther downstream the helical nature of the large-scale 
structure is dominant (Viswanathan and Sankar 19 ). Also, 
both experiments and hydrodynamic stability analysis 
suggest that the m = +1 helical mode is dominant above a 
jet Mach number of about 1.4. Thus by including 
azimuthal mode (m = +l) instability in a three- 
dimensional analysis, predictions of sound pressure levels 
could be considerably improved (Garnet and 
Estivalezes 20 ). Such three-dimensional simulations require 
considerably more computing resources, which are 
outside the scope of the present work. Further 
improvements in the predictions can be realized by 
considering inflow disturbances of various frequencies (or 
random disturbances). 

Figure 1 1 shows the corresponding comparison with 
the OASPL plotted as a function of the axial distance 
from the exit plane. The computed results of OASPL for 
three values of z/r ; = 30, 55 and 80 are displayed in 

figure 12, the last one corresponding to the data by 
Seiner 14 . With increasing values of zlrj , the peak value 

of OASPL decreases, and the corresponding peak angle 
decreases. This is to be expected in view of the 
directivity of Mach wave radiation. The correctness of 
this trend also serves to verify the Kirchhoff code 
implementation. 

The spectral content of SPL (dB/Hz) along a 
direction corresponding to the maximum radiation 
(x/rj =60 ,z/rj =80) is displayed in figure 13. Figure 

14 indicates the variation of 1/3 octave SPL (dB) as a 
function of the frequency. A well-defined peak SPL value 
of about 128 dB is seen to occur at a frequency of about 2 
kilohertz (kHz), which corresponds to a Strouhal number 
of 0.2. This value is consistent with the frequency of the 
imposed disturbance at the nozzle exit. Figure 15 
presents the result of figure 14 in a semi-log plot. 

6. Conclusions 


The variation of the OASPL in the streamwise 
direction at z/r ; =80is compared in figure 10 with the 


The CFD/Kirchhoff approach is capable of predicting 
the SPL spectrum and the OASPL at any location in the 
near-field and the far-field. The relative reduction in the 
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computational time due to the Kirchhoff method and the 
axisymmetric simulation make it an attractive engineering 
tool. The one-equation turbulence model of Spalart- 
Allmaras for computing turbulent stresses affords 
additional reduction in the computational time and at the 
same time improves the robustness of the code relative to 
two-equation turbulence models. The peak angle of Mach 
wave radiation and the peak level of OASPL are 
satisfactorily predicted. Predictions at off-peak 
conditions could be substantially improved by considering 
three-dimensional disturbances (azimuthal mode 
instability) and disturbances of various frequencies. 
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Table 1. Summary of Nozzle Conditions 


Stagnation temperature, K 

755 

Nozzle exit diameter, m 

0.09144 

Exit pressure, atm 

1.0 

Exit temperature, K 

419 

Exit velocity, m/s 

822 

Acoustic velocity at exit, m/s 

41 1 

Nozzle exit Mach number 

2.0 

Exit jet Reynolds number 

1.3x10° 

Ambient pressure, atm 

1.0 

Ambient temperature, K 

280 


z = 3.66 m (Langley 



Figure 1. Schematic of the Test Case Nozzle 
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Figure 2. Axisymmetric grid (200 x 100) 



Figure 3. Mean Axial Velocity Contours in the Jet From CFD Solution 
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Figure 4. Jet Centerline Variation of Axial Velocity From CFD 
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Jet Centerline Variation of Static Pressure From CFD 
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Figure 6. Jet Centerline Variation of Mach Number From CFD 
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Figure 7. Instantaneous Pressure Contours in the Jet From CFD Solution 
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Figure 8. Contours of OASPL in the Source-Field as Computed Directly From CFD 
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Figure 9. Contours of OASPL Predicted From CFD/Kirchhoff Method 
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Figure 10. Comparison of Angular OASPL Distribution at z! r } =80 
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Figure 11. Comparison of Axial OASPL Distribution at z I =80 



Figure 12. Computed Axial Distributions of OASPL for z / r- =30, 55, and 80 
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Figure 13. Computed SPL Distribution at x! r. = 60, and z/ r. = 80 
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Figure 14. Computed SPL Distribution at x! r ; - - 60, and z/ r, = 80 
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Figure 15. Computed SPL Distribution (Semi-Log Plot) at x! r. 


60, and z! r ; = 80 


14 


